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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

ANALYSIS AND PRELIMINARY INVESTIGATION OF EDDY-CURRENT HEATING 
FOR ICING PROTECTION OF AXIAL-FLOW-COMPRESSOR BLADES 
By Thomas Dallas and C. Ellisman 


SUMMARY 

An analysis and a preliminary experimental investigation of 
the use of eddy currents for heating the blades of an axial-flow 
compressor of a turbojet engine as a means of ice protection is 
presented. Steinmetz's analysis of eddy currents in flat plates 
is extended and formulas are derived that permit calculations of 
currents, voltages, flux distribution, and power when operation 
is not limited to conditions of constant permeability. 

The application of eddy-current heating appears to be fea- 
sible for the generation of power densities required for icing 
protection of axial-flow-compressor blades. 


INTRODUCTION 

The advent of all-weather flying of jet-propelled aircraft 
demands provision for complete ice protection of the engines. 

One component of this type of engine that has not been provided 
with ice protection is the compressor inlet, although NACA flight 
investigations of an axial-flow-compressor-type turbojet engine 
have indicated that the compressor inlet is subject to icing. 

Three methods have been proposed by the NACA to protect the 
engine: (l) Inertia separation of the water droplets from the 

induction air (reference 1); (2) heating of the charge air 
(reference 2); and (3) heating of the engine parts by hot gas 
(by a method similar to that of reference 3). A fourth method 
of de-icing, which involves heating of the engine parts by the 
use of electrical eddy currents that can be induced if a pul- 
sating magnetic field is created in the blades to be protected, 
is proposed by the NACA Lewis laboratory and discussed herein. 

An analysis of eddy currents in a compressor blade, the 
geometry of which has been limited to that of a flat plate for 
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simplicity, is presented. Formulas are derived that permit cal- 
culations of currents, voltages, flux distribution, and power when 
operation is not limited to conditions of constant permeability. 
Preliminary investigations are described in which a model blade is 
subjected to flux variations and the heat generated is determined 
as a function of magnetizing ampere turns. The heat required for 
adequate protection of the inlet guide vanes is estimated. Con- 
clusions based on experiment and theory are stated regarding the 
suitability of the method, and statements are included regarding 
the significant parameters. 


DESCRIPTION OF EDDY-CURRENT HEATING 

The configuration of a typical axial-flow compressor (fig. l) 
lends itself to a design that incorporates an eddy-current gen- 
erator. The principle of eddy-current heating as applied to the 
inlet guide vanes is illustrated by the schematic diagram of fig- 
ure 2. In order to apply eddy-current heating, the blading involved 
must be of ferromagnetic material. In addition, a ferromagnetic 
flux path must be provided to complete the magnetic circuit. A 
direct-current coil is needed to supply magnetomotive force required 
for operation of the system. A flux-chopping mechanism (section A-A, 
fig. 2) provides a varying air gap. The varying reluctance of the 
air gap causes a flux change, which results in the creation of eddy 
currents and heat. 

It will subsequently be shown that it is desirable to use flux 
variations of high frequency. In a typical axial-flow compressor 
having 88 inlet guide vanes and operating at 7600 rpm, a cursory 
examination (without investigation of all design criterions) indicates 
that 44 chopper teeth could be used to give a frequency of 5600 cycles 
per second. At this frequency, any appreciable flux variation in 
parts of the magnetic circuit other than the blading would cause 
undesired heating and loss of energy. In addition, if flux variation 
in parts other than the blading is allowed, the time constants of 
the coil and other parts involved may make it impossible to obtain 
the full amplitude of flux variation in the blades at the frequencies 
under consideration and the system may fail to operate. (Because 
it is impossible to maintain constant flux in all parts of the chop- 
per, lamination of the chopper will minimize the difficulties 
discussed. ) 

The foregoing principles dictate a design criterion typically 
illustrated in section A-A of figure 2. The teeth in the chopper 
should be of such a length that a6 one blade is leaving a tooth 
another blade is taking its place in the magnetic circuit. In this 
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manner, the only parts of the magnetic circuit that are experiencing 
a flux change are the blading and the immediately adjacent parts. 

A constant flux is thus maintained in the rest of the magnetic 
circuit and a constant flux always links the coil. 

In a general application of eddy-current heating where more 
than one set of blades is protected, the chopper may be either in 
the rotor or the stator but the blades to be heated must pass over 
the teeth of the chopper. A chopper must therefore be provided 
for each row of blades to be protected. 

In compressor -blade heating by means of eddy currents, flux 
in the blades can be either unidirectional-pulsating or alternating. 
In the following analysis, however, only the alternating- flux case 
is considered. An experimental correlation between this type of 
operation and the unidirectional-pulsating type is presented. 


SYMBOLS 


A 

®max 

®max 


c 

d 

E 

®max 

f 

•^max 

I 


The following symbols are used in this report: 

cross-sectional area of blade, square centimeters 
maximum value of alternating flux density, gausses 

u maximum value of alternating flux density under uniform- 
flux distribution, gausses 

constant, /\/o.4 f ^ io~®, centimeters - ^ - 

plate thickness, centimeters 

root-mean- square value corresponding to E maT , volts per 
centimeter 

maximum value of alternating voltage gradient (in 
direction of plate width), volts per centimeter 

frequency, cycles per second 

maximum value of alternating-magnetic-field intensity 
(in direction of plate length), oersteds 

root-mean- square value corresponding to X majc , amperes 
per centimeter 
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^max 

maximum value of alternating current per centimeter length 
of plate; equilibrant of all currents in incremental 
vo limes included from plate center to any plane parallel 
to outer plate surface, amperes per centimeter 

j 

complex operator equivalent to A /^l 

L 

length of blade, centimeters 

l 

distance from plane ABC (fig. 3) toward outer surface of 
plate, centimeters 

l o 

one-half of plate thickness d, centimeters 

l v 

depth of penetration at constant permeability |i, 
centimeters 

*p av 

depth of penetration at average permeability n av , 
centimeters 

z T 

total depth of penetration + Ij, centimeters 

*2 

upper limit of Z for given impressed field intensity 
of Bmax 2> centimeters 

NI 

total root-mean- square ampere turns 

^a^max a 

maximum value of alternating ampere turns required across 
air gap 

^b^max b 

maximum value of alternating ampere turns required across 
blade 

P T 

total air-gap permeance, centimeters 


Pl» . .Pg permeances of various volumes composing air gaps, 
centimeters 


P 

eddy-current power generated for uniform flux distribu- 
tion, watts 

p e 

eddy-current power generated per square centimeter of 
blade surface, watts per square centimeter 

p e' 

eddy-current power generated per square inch of blade 
surface, watts per square inch 


1125' 


1125 


NACA EM E9E06 


5 


0 o 

n s 

distance from surface ABC (fig. 3) toward center of plate, 
centimeters 

V 

plate volume, cubic centimeters 

a 

mean of 3-& (fig* 6), degrees 


limiting value of a, mean of 3 2 -& 2 > degrees 

3 

angle between I mBy ^ and l max (fig* 6), degrees 

P 2 

upper limit of 3, angle between I^y 1 and I mnY 
degrees 

7 

conductivity, mhos per centimeter 

6 

angle between 3> mn y x and <l' max (fig. 6), degrees 

5 2 

upper limit of 6, angle between $ . and $ 0 , 

degrees 

e 

power-factor angle (fig. 5), degrees 

n 

constant magnetic permeability, (B mny . i /B^y 

^av 

average magnetic permeability 

P 

resistivity, (l /y) , ohm-centimeters 

$> 

max 

maximum value of alternating flux included in cross- 
sectional area having width of 1 centimeter and thick- 
ness with limits from plate center line to an outer 
limit, which is plane where field intensity is H^y, 
maxwells 


Subscripts: 


s 

values in or to plane ABC of plate (fig. 3) 

1 

values in or to plane ABC (fig. 3) when field intensity 
at plane ABC has reached value at saturation point 
H max ^ on magnetization curve 

* 2 

upper limiting values where field intensity is given 
impressed field intensity H^y 2 and flux density 

is saturation value B TO „„ 

max 1 


Overscoring indicates vector quantities. 
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ANALYSIS 

Development of Analysis 

In eddy-current calculations for plates, the following formula 
is commonly used for determining power when a sine wave of alter- 
nating flux at frequency f is employed: 

10 -16 watts (1) 


This expression is valid if the magnetic effect of eddy currents is 
negligible and if the lines of force are uniformly distributed 
throughout the plate and are parallel to the sides of the plate. 

Eddy currents, by virtue of their countermagnetomotive force, 
tend to damp the flux in the center of the cross section of the 
magnetic material. The depth of penetration of the flux into the 
surface of a plate of magnetic material is given by Steinmetz 
(reference 4, p. 371) as 



3570 

pav= /s/w 


( 2 ) 


At a frequency of 6100 cycles per second, for instance, an 
assumed average permeability of 119, and a conductivity of 10'" mho- 
centimeters corresponding to soft iron, the penetration would be 
0.0132 centimeter or 0.0052 inch. With an average blade thickness 
of approximately 0.062 inch, it is apparent that because of the 
damping action, the flux density is negligible throughout the 
greatest portion of the blade cross section and that a serious 
error is involved in considering it otherwise. 

Steinmetz (reference 4) has considered alternating-magnetic- 
flux distribution only under conditions of constant permeability. 

In the present work, however, operation well past the knee of the 
magnetization curve is necessary and the work of Steinmetz is not 
directly applicable. Steinmetz considered a section of a flat 
plate or lamination similar to that shown in figure 3. He divided 
the plate into incremental slices of width ds at distance s 
from the surface of the plate and set up the partial differential 
equations of the system. The direction of the flux was assumed 
along the length or longest dimension of the plate. 
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Steinmetz found that, if the flux at_the center of the plate 
is negligible, the maximum flux density Bn, ST at any point at 
depth s from the plate surface is given by 

*W = B max s e' CS (cos cs - J sin cs) (3) 

An examination of equation (3) reveals the nature of flux 
penetration into the surface of the plate. At the surface of the 
plate, or when s = 0, is equal to B^y g , which is the 

flux density produced by the impressed magnetomotive force. The 
real part of equation (3) is therefore in time phase with the 
impressed magnetomotive force. It is apparent that the flux lags 
in time and decreases rapidly in value with depth s into the 
plate. At a depth equivalent to cs equal to jt/2, the flux 
lags 90° and is decreased to 20 percent of its value at the sur- 
face. When cs is equal to «, the flux lags 180° but its 
intensity is only 5 percent. At the point where the flux is again 
in phase with the impressed magnetomotive force, its value is only 
0.2 percent. 

Integration of equation (3) between the limits s = 0 and 
s = Z Q produces the maximum flux in one-half of the plate for 
1 centimeter of width. The condition that the flux in the center 
of the plate is zero dictates that e -c ^o is negligible; therefore 


max s 


B, 


max s 
2c 


(1-J) 


(4) 


or 


$ 


max s 


W s 
Ji c 


(4a) 


Equation (4) may be used to obtain the mean maximum flux density 
Bmax m b 7 dividing by Z Q : 



^max s 
2ci 0 


(1-j) 


(5) 


If m is made equal to B ltlfl7 g , then 2 Q becomes the depth 

of penetration Zp or is the depth to which the flux would pene- 
trate to produce the same total flux if the flux density at the 
surface was maintained throughout the section of the plate. 


8 


NACA EM E9E06 


1 p - is 


( 6 ) 


or 


l 


1 


( 7 ) 


P ^ c 


Inasmuch as Steinmetz's work is not directly applicable because 
of operation past the knee of the magnetization curve, further 
development is necessary. Consideration will still be limited to 
the flat plate of figure 3. No generality will be lost by this 
limitation on shape, however, because effects of eddy currents will 
be shown to be independent of plate thickness d and to be func- 
tions of the surface only (provided flux density at the center of 
the object is negligible). 

A brief general statement will establish the concept of the 
analysis. Until saturation is reached, the absolute value of the 
maximum flux density B max s in the outer surface increases with 
impressed magnetomotive force and decreases in absolute value 
exponentially with distance into the material according to equa- 
tion (3). As the impressed magnetomotive force is increased beyond 
saturation, the maximum flux density in the outer surface cannot 
increase further, but a different action results. The absolute 
value of the flux density remains essentially constant for a short 
distance into the plate before decreasing exponentially. The dis- 
tance of this penetration at constant flux density is determined 
by the amount of impressed magnetomotive force in excess of that 
required for saturation. 

The normal magnetization curve for SAE 1020 steel, annealed at 
1500° F and furnace- cooled, was chosen as a representative curve 
for purposes of analysis. The curve, which is plotted in figure 4 
from data of reference 5, is a plot of flux density in gausses as 
a function of field intensities in oersteds. The magnetization 
cui*ve is approximated by two straight lines; the line AF drawn from 
the origin approximately tangent to the curve and another line 
representing the saturation flux density B mav -|_ drawn parallel to 

the abscissa axis at approximately the mean flux density for a 
given range of field intensity. The mean is estimated from the 
geometry of the curve of figure 4. As will subsequently be shown, 
most of the heat is generated in the outer layers of the plate at 
the higher magnetomotive forces, which is justification for giving 
greater weight to the higher flux densities. The choice of the 
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mean is not critical because i enters into the final power 

formula to the one-half power and a 4-percent error in estimating 
Bmay ]_ results in a 2-percent error, approximately, in the power. 
The lines BC, BE, and FG are for limiting values of field intensity 
of 6, 30, and 50 oersteds, respectively. In order to make the 
problem specific, let the limiting maximum field intensity H^y 2 
be 30 oersteds; then the curve will be represented by lines AB and 
BE. It is apparent that Steinmetz's equations are applicable along 
the line AB where the permeability q, which is the slope of 
line AB, is constant. 

Equations (3) to (7) can be represented by the vector diagram 
of figure 5. Time phase is measured by counterclockwise rotation. 
The reference vector I^y s extends from the pole toward the 
upper left. The equality of the flux components of equation (4) 
fixes the current at a lagging power-factor angle of 45° with 
respect to the impressed voltage. The voltage of self-induction 
E max lags its flux g by 90°; whereas the impressed or 

consumed voltage, which is equal but opposite to E rnftV s , leads 

by 90°. This diagram can be chosen to represent conditions at the 
saturation point B of figure 4. Then other points along the 
line AB, which represent conditions at deeper penetration into the 
plate, can be represented by vectors having the same relative pro- 
portionality and phase, but which are shifted clockwise to the 
proper lagging time phase relative to conditions at point B and 
which are reduced in magnitude by the factor e" cs . 

Two more equations are required to establish the magnitudes 
of the voltages and current in figure 5. The maximum voltage 

of self-induction per centimeter width of plate from the 


In order to eliminate any turns-ratio numbers, the assumption is 
made that the impressed magnetomotive force is supplied by a one- 
turn coil. As previously stated, the flux density B^y g at the 

outer surface of the plate is that due to the primary magnetomotive 
force alone because no effects of bucking eddy currents exist at 
the surface. Because the flux is assumed to be along straight 
lines parallel to the sides of the plate, the equation for the 
intensity of the field at the surface is 



max s 


is 


Emax s - $max s ^ 


-8 


( 8 ) 


^max s ~ n -^max s 


(9) 
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where 

Ifflpy g value in SYMBOLS (also, current required in a one- turn 

energizing coil to furnish a field intensity of B mnx g 

for volume of plate (fig. 3) 1 cm wide, 1 cm long, and 
depth equal to one-half plate thickness), amperes per 
centimeter 

The vector diagram of figure 5 represents the flux, the current, 
and the voltage associated with the volume of the plate described 
in the preceding symbol definition. With the vector diagram of 
figure 5 chosen to represent the conditions at the point D of the 
magnetization curve (fig. 4), the outer surfaces ABC and A'B'C' of 
figure 3 will be saturated at a flux density ^ ^ under a field 
intensity of B^y Equations (4a) and (7) now become 


<J> 

max 


1 


^max 1 
ife c 


(4b) 


and 



(7a) 


In the analysis it was found convenient to use the following 
physical concept. The assumption was made that thin sheets of 
SAE 1020 steel of incremental thickness dZ are added to the sur- 
faces ABC and A'B'C' and become integral with the original plate 
of material. In addition, it was assumed that the impressed mag- 
netomotive force is increased to balance the bucking eddy currents 
in the new sheets in order that the flux densities at planes ABC 
and A'B'C' remain unchanged. Then the vector diagram of figure 6 
will represent the new conditions. 


The dashed vectors of figure 6 are a duplication of figure 5, 
with the exceptions that the flux- component vectors have been 
omitted for clarity and the vectors now represent conditions at 
the saturation point. The diagram continues the representation 
of a plate volume having a surface 1 centimeter square and a depth 
equal to one-half the plate thickness. Regardless of saturation, 
it is assumed that currents, voltages, and flux are still repre- 
sentable by vectors and, in the case of figure 6, by maximum 
values of sine waves. As the new sheets of thickness dZ at 
constant flux density B mflT ^ are added, the vectors will change 
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in absolute value and position, as indicated by the dashed lines 
representing their loci. It was assumed that sufficient sheets 
have been added to reach the limiting maximum field intensity 

g 30 oersteds that was chosen. The various quantities at 

this condition will be identified by the subscript 2. 

When a new sheet is added, the total flux is increased by the 
vector q di, which is in phase with the current I ^y , as 

shown in the diagram (fig. 6). From equation (4b) and with B^y 

a constant, it is possible to write the equation for the absolute 
value of the maximum flux from the diagram 

«W = Vz 1 * / "W 1 c ° s < 45 ° . 5«a + 

Jo 


Bn, n y x Z cos (45° +a) (10) 

where a is the mean value of (3-6) for the given value of Z. 

The infinitesimal change in maximum field intensity d H rT1fly 

between the sides of a sheet is proportional to the current in the 
sheet 


or 


d Vx = °* 4 * ? W dl 


^ax 

dZ 


= °* 4 * 7 Emax 


( 11 ) 


where E [[iay is the maximum voltage per centimeter width induced in 
the sheet by the included flux $ m „y The relation between E^y 


and $ 


max 


is 


Ejnax = 2 jr f $ 


max 


10 


-8 


( 12 ) 


From equations (10), (ll), and (12), the following relation is 
obtained: 
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fflmax et 0A * Z Vi 1 

dZ 4 c 


& + 


cl cos (45° + a.) 10 


a)] 


-8 


(13) 


When ^ \ny i = B max 1 substituted an< l equation (13) is inte- 
grated in the limits from 0 to Z 2 , to 2 > the 

absolute value of the maximum field intensity 2 is 


\ax 2 = B max 1 + c H max 1 


cl. 


1 2 + a/2 — g — cos (45° +a^) 


(14) 


from which the value of Zg will be 


z 2 = - 


1 + 1 1 


^max 


^2 c cos(45°+a 2 ) 'Nj 2c 1 2 cos 2 (45°+a 2 ) c 2 cos (45° +a 2 )y H max i 

(15) 


or 


1 + t/2 cZ 2 cos (45° + 02 )= aJ 1 + 2 cos ( 45 ° +a >) - 1 


Also, squaring both sides. 


fj2 cZ 2 + c 2 Z 2 2 cos (45° + a 2 ) = ~ ~ ~ c 

nnax 1 


^nax 2 


- 1 


(15a) 


(15b) 


The total depth of penetration is the sum of equations (15) and (7a) 


1 T = I P 


l +l 2 


cj2 


+ 


cos(45 c 



1 + 



(15c) 

The voltage gradient E mflT 2 in the outer plate surface can be 

derived by using equations (10) and (12) and substituting the upper 

limits $ for $ , and l 0 for Z, 

max 2 max & 
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2 = a/ 2 * f 10" 8 Em ^ - ~ |l + s[2 cl 2 cos (45° -Wg)] (16) 

Substitution of equation (15a) and division by a/ 2 to obtain the 
effective value Eg on the assumption of a sine wave, gives 

Eg • it f 10 -8 - m ^ - 1 [\jl+2 cos (45° + a 2 )Qj^^ - 1 ^ (16a) 


where 

E 2 value in SYMBOLS (also the root-mean- square value of impressed 
voltage required for plate volume having surface 1 cm square 
and depth equal to one-half plate thickness), volts per 
centimeter 

The relation of equation (15a) may now be used to eliminate l 
from the equation for the maximum flux Equation (10) may be 

rewritten by going to the limits 2 for and 1 2 for Z 

g 

$ max 2 “ [} + ^ cl 2 cos (45°+a^j (10a) 

or, when equation (15a) is used, 

*« * - ^ 2 cos( 45 °^ ife' 1 ) (10k) 

The argument involved for obtaining the current 1^ 2 is 
identical to that used for equation (9); the replacement therefore 
of Lt s B max 2 in e< l ua ' bion ( 9 ) gives the expression for 

I max 2 

W 2 - °- 4 * W 2 (17) 

where 

2 value in SYMBOLS (also, current in amperes required in one- 
turn energizing coil to furnish given impressed field 
intensity 2 for volume of plate having 

surface 1 cm square and depth equal to one-half plate 
thickness), amperes per centimeter 
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is 


The effective value of current Ig for an assumed sine wave 


h = 


W 2 _ ^max 2 
^2 V2 0.4jc 


(17a) 


The remaining derivation is that of the power-factor angle, 
which will be obtained by finding first the angle 6g and then Pg. 

The derivation of 8 2 and Pg may be found in appendix A. 


B 2 = 


tan (45° +a^> 


log c 


and 


Pg = cot (45° +a 2 ) log e 


1 ♦ 2 cos - l) 


1 . sin (45°« 2 )fe - l) 


(18) 


(19) 


For a given value of H maz 2 /^nax 1 > various values of ctg 

were substituted in equations (18) and (19) until, by successive 
approximations, an average value of cig was obtained that satis- 
fied the equation 


“2 



(P 2 -6 2 ) d(p 2 ) 


where (p«) is the limit of integration. The resulting values 

Cm q, 

of P 2 ~8 2 are shown as a function of P 2 in figure 7. The values 

of -®max 2 /^nax 1 an< ^ a 2 are P os ^ e< ^ heside each calculated 

point. Eeplots of the data from figure 7 are presented in fig- 
ures 8 and 9. In figure 8 the angle ctg is presented as a func- 
tion of various values of H ^y 2 /Hmax and figure 9 shows the 
angie Pg-6g as a function of 2 / H max 1* 


Results of Analysis 


The eddy-current power p e absorbed as heat per square centimeter of plate surface 
will be the product of the voltage, the current, and the power factor (fig. 6) and may be 
expressed as 

p e = E 2 I 2 sin (45°+|3 2 -5 2 ) 

or from equations (16a) and (17a) 


Pe . f(10 -S) ^ ^ .1. 


Revision yields 

10 ' e 


r& 


^e " 


1+2 cos (45°-HXg) ^— — - ij sin (45°+8 2 -6 2 ) 


(20) 


Replacement of c by its parts yields 


Pe=' 


1 1 

1.25 2_2 


f- 5 ^ Vl2Ul 10- 4 A All-2008 (4tf>«2j| +2 008(450402)5!!^ sio (4S°4(S 2 -S 2 ) 

V max 1 


watts per square centimeter (20a) 


where p is resistivity (ohm-cm). Multiplying equation (20) by 6.45 gives the equation 
for the heat generated per square inch. 


11.4 f 


v- 


®max 


2 ®max 1 


10 


-8 


/l+2cos ^5 0 +a2> 


max 2 

a 

max 1 


- 1) sin (45°+3 2 -8 2 ), watts per square 


(20b) 


En 


inch 
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Operation at Constant Permeability 

The equations that have been derived are not limited in use to 
conditions in which the outer plate surface is saturated, but are 
also applicable in the region of the approximated magnetization 
curve where the permeability p. is constant. At the saturation 
point, where 2 is e< l ual to Ejnax 1' tlie e ! uation for the 

total depth of penetration (equation (15c)) reduces to equation (7a). 
The depth of penetration remains constant at the value given by 
equation (7a), even if the impressed field is reduced below ^ 

The power-factor angle 9 remains 45° for field intensities equal 
to or less than H^y and p 2 and vanish. At the satura- 

tion point, the radical in the voltage equation (16a), the total- 
flux equation (10b), and the power equation (20b) reduces to a 
value of 1. These equations are still applicable for impressed 
field intensities less than H^y with the substitution of the 

impressed field intensity g for H max 2 and the value 

equivalent to the impressed field s for B mn y ^ after the 

radical has been replaced by 1. The total-flux equation (10b) 
reduces to equation (4a) . The current formula (17a) applies if the 
given impressed field intensity H max is used for H^y 2 . 


EXPERIMENTAL INVESTIGATION 

Prior to any attempt at full-scale application of eddy-current 
heating of the inlet guide vanes, preliminary investigations were 
conducted from which information was obtained as to the number of 
field ampere turns required to produce specified amounts of heat. 

The schematic layout of the equipment in which sample blades 
were investigated is shown in figure 10. For simplicity, instead 
of using a chopper to vary the flux, a similar effect was accom- 
plished by varying the magnetomotive force or ampere turns by con- 
necting the coil shown in figure 10 to a 500-watt, 6100-cycle-per- 
second, audio-frequency power generator. In order to match the 
load to the generator properly, the turns on the coil were varied 
from 316 to 210, in addition to resonating this load with series 
capacitors. The current in the coil was measured with radio- 
frequency ammeters. 

The heat generated was determined by measuring the tempera- 
ture rise and the rate of water flowing through the jacket. Heat 
exchange to the laminated core was minimized by keeping the core 
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temperature near the mean temperature of the water. From runs with 
and without insulation, insulation on the water jacket to reduce 
heat transfer between the atmosphere and the jacket apparently was 
unnecessary in this investigation. 

The pulsating unidirectional operation was synthesized by 
superimposing a direct current and an alternating current in the 
inductor of the experimental unit, as shown in figure 11. The 
direct current to the inductor was introduced by means of a filter- 
ing network. The 0.55-henry choke, which was especially wound with 
a large air gap to prevent saturation, had approximately 10 times 
the inductance of the inductor and thus shunted only about 10 per- 
cent of the alternating current from the inductor. The inductance 
value of the remaining choke was not critical, but its impedance 
was sufficient, relative to the 20-microfarad condenser, to force 
most of the alternating current of the 0.55-henry choke through 
the measuring meter in series with the 20-microfarad condenser. 
Because the ratio of inductive reactance to resistance of the 
0.55-henry choke and the inductor are sufficiently alike to make 
the currents in phase for all practical purposes, the difference 
between the readings of the two alternating- current ammeters was 
therefore the alternating current in the inductor. No appreciable 
alternating current passed through the direct- cur rent ammeter or 
the battery circuit and the direct- current meter read the true 
direct current applied to the inductor. 

The first investigations were made on annealed SAE 1020 steel 
with a total air gap between the laminated iron core and a sample 
blade of approximately 0.04 inch. The air gaps used for these 
investigations were approximately those that would exist in prac- 
tice between the blade and the adjacent part through which the 
flux would flow. In order to permit a correlation between experi- 
mental results and the analysis, the blades were made of rectangu- 
lar cross section. The actual size of the sample, however, which 

was 0.125 by 1 by 4^ inches, is approximately the same as the inlet 

guide vanes of a typical axial-flow compressor. The second and 
third experiments were made with a blade of Armco Magnetic Ingot 
Iron of the same size. In the second experiment, the air gap was 
kept at approximately 0.04 inch and in the third experiment the 
air gap was increased to approximately 0.06 inch. A sample of 

Hipernik, 0.094 by 1 by 4| inches, with a total air gap of 

0.04 inch was used in the fourth experiment. This material has 
high resistivity and high maximum flux density. The investiga- 
tions on all samples were made by varying the number of magnetiz- 
ing ampere turns and measuring the heat generated. 


18 


NACA EM E9E06 


The fifth experiment was made with the synthesized pulsating 
unidirectional current. As illustrated by the curve in figure 11, 
the two currents were apportioned to subject the blade to a magneto- 
motive force that varied from 0 to 2 a/2 times the root-mean- square 
value of alternating magnetomotive force at each experimental con- 
dition. The Armco iron sample was used in this investigation and 
the total air-gap spacing was 0.04 inch. 

A 10-turn search coil was placed around the center of the 
sample blade in all experiments for the purpose of measuring the 
voltage and viewing the voltage wave shape on a cathode-ray 
oscilloscope. Ten turns were used to obtain sufficient signal 
voltage to make the hum voltage and other spurious voltages 
negligible in comparison. 


CALCULATION OF AMPERE TURNS 

The correlation between the analytical equations and the 
experimental results depends on several new calculations that must 
be made in addition to computations using the equations already 
derived. 

The blade was positioned in the experimental rig (fig. 10) 
to make the air gaps at each end of the blade approximately equal. 
Figure 12 shows the geometry used in calculating the permeances at 
each blade end. No attempt was made to draw figure 12 to scale or 
to make the relative sizes of the various volumes correct. (Volumes 
are designated by numbers 1 to 8.) The purpose of the figure is to 
give the geometric shapes of the various volumes to illustrate the 
application of formulas from reference 6. The suggestion given in 
this reference was followed in limiting the fringing-flux calcula- 
tions to the entire space immediately surrounding the pole faces 
and not at a great distance along the blade (h+t+g was limited 
to approximately 0.75 in., where h, t, and g are dimensions, 
as shown in fig. 12). The geometry of figure 12 is based on the 
two propositions that the flux lines are normal to the metal sur- 
faces and that such surfaces are equipotential . These propositions 
are obviously assumptions, because the magnetomotive-force drop 
along 0.75 inch of blade is appreciable. 

The formulas for the permeances of the various volumes desig- 
nated by the numbers in figure 12 follow for convenience. The air 
gap directly under the blade was assumed to have no flux, except 
for a volume equivalent to the depth of penetration. This per- 
meance was therefore taken to be equal to the product of the 

depth of penetration Irp and the circumference of the blade cross 
section, divided by the distance g: 
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22 T (d+w) 

5 


where w refers to width, as shown in figure 12. Also 

P 2 = 0.52 w 



The plane ABC is so positioned that the average value of volume 5 
remains unchanged. In this manner, volumes 4 and 5 can he com- 
bined in one calculation. On a scale drawing, the line BA was 
found to make an angle of 50° with the horizontal, the angle ABD 
therefore being 140°. Then 

•o / 180\ w n ( g+h+t\ 1.28 w„__ (e+h+t\ 

P 4 + 5 ” (jtoj H log e V 8*7 — log e \^r) 

and 


P 6 = 0.52 d 



Volume 8 is considered a partial quadrant of a spherical shell of 
thickness ..t+h. The fraction used as a factor in making the com- 
putation is the ratio between 180 and the arithmetic mean of 140 
and 90, or 


•n 180 t+h « iq / j . 

P 8 ~ ^ 140+90 j 4 “ °* 39 ( t+h ) 


The total permeance for one end will be 


2 (-T + P 2 + ?3 + P 4+5 + P 6 + ^7 + ^8 
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The total air-gap permeance will equal one-half of this value, or 


P T = ~z + P 2 +P 3 +P 4+5 +P 6 +P 7 +2P 8 


( 21 ) 


The following equation can be written from the total- air- gap 
magnetomotive -force equation for the real component 


N I 

a max a 


0.4 it l 


A ^max 2 cos (45°+P 2 -& 2 ) 


Substitution of equation (10b) gives the equation for air-gap 
ampere turns 


N I 

a max a 


cos(45°^ 2 -& 2 )B max 1 


0.4 nhJZ cl , 


•T 




( 22 ) 


The final determination needed is that for the magnetomotive- 
force drop along the blade. Previous discussion showed that the 
flux density at the outer surface of the blade is the result of the 
impressed magnetomotive force alone; the following equation is 
therefore true: 


Thus 


_ ft -^b^max b 
■“max 2 L 


N. I 

b max b 


L H . 
max 2 

0.4 « 


(23) 


The total effective ampere turns required are 


NI 


N I 

a max a 


K. I , 
b max b 


V 2 


(24) 


When the previously derived equations are used, the eddy- 
current power per square inch of blade surface appearing as heat- 
can be calculated from equation (20b), and the voltage can be 
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determined from equation (16a). The total depth of penetration 
of the flux into the blade is found from equation (15c). 

Values of ^ 2 selected for making check calculations were 

6, 30, and 50 oersteds. These values are shown in figure 4 for 
SAE 1020 steel, together with the corresponding values of 11,500, 
15,000, and 15,800 gausses for the average values of -|_. The 

permeability of the left part of the magnetization curve, which is 
the slope of the line AF, is 2665. The magnetization curve for 
Armco Magnetic Ingot Iron (reference 7) and lines BC, DE, and FG 
indicating the average values of ^ chosen for 6, 30, and 

50 oersteds, respectively, are shown in figure 13. 

The magnetization curve for Hipernik (from reference 8), which 
has relatively high resistivity, is shown in figure 14. The hori- 
zontal lines BC, DE, and FG on this curve represent chosen average 
values of B ma;j!; 1 of 12,800, 14,000, and 15,000 gausses for limiting 

field intensities 2 of 6, 20, and 40 oersteds, respectively. 

RESULTS AND DISCUSSION 
Electrical Considerations 

Results of the calculations and the experiments are given in 
table I and plotted in figures 15 and 16. Figure 15 is a plot of 
voltage at the surface per centimeter width of blade against effec- 
tive impressed ampere turns. The search coil around the center of 
the blade is the source of the voltage data, which are divided by 
the product of coil turns and blade cross-sectional circumference 
in centimeters to obtain the values for the plotted points. The 
curves are for investigations on SAE 1020 steel, Armco Iron, and 
Hipernik, each with total air gaps of 0.04 inch. Calculated values 
of voltages from equation (16a) for field intensities at the blade 
surface of 6, 30, and 50 oersteds are also plotted on figure 15 for 
Armco Iron and SAE 1020 steel. The limiting field intensities used 
in making voltage calculations for Hipernik were 6, 20, and 40 oer- 
steds. The root -mean- square or effective ampere turns for each inten- 
sity were calculated from equations (24), (23), and (22). The purpose 
of these data is to show the close agreement between calculated and 
measured voltages. The numerical comparison between calculated and 
measured voltages obtained from columns 10 and 11, respectively, of 
table I shows agreement to within 10 percent, except for Hipernik 
where the disagreement is approximately 16 percent. The magnetiza- 
tion data used in the calculations, however, came from published 
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average curves that may not exactly represent the samples inasmuch 
as the measured conductivity of the samples used in the calcula- 
tions did not agree with the published average data. 

Three of the curves in figure 16, which represents eddy-current 
power generated per square inch of blade surface plotted against 
effective impressed ampere turns, are for SAE 1020 steel, Armco 
Magnetic Ingot Iron, and Hipernik for a total air gap of 0.04 inch. 
Figure 16 shows the close correlation between the curves and plotted 
check points obtained from analytical formula (20b) for field inten- 
sities at the blade surface of 6, 30, and 50 oersteds for Armco 
Iron and SAE 1020 steel and 6, 20, and 40 oersteds for Hipernik. 
Again the root -mean- square ampere turns for each of these inten- 
sities were calculated from equations (24), (23), and (22). From 
columns 12 and 13 of table I, the agreement between calculated and 
measured values is within 10 percent. 

On the Armco Iron curve for a total air gap of 0.04 inch 
(fig. 16), the data of the synthesized-pulsating-flux investiga- 
tions are plotted. The measured eddy-current power is plotted 
against the root-mean- square value of the alternating component of 
current. These points agree closely with the simple alternating- 
current tests and indicate that the superposition of the direct- 
current magnetomotive-force component has very little effect if 
any. It is important to understand, however, as the curve in fig- 
ure 11 demonstrates, that for any given power the ampere turns from 
figure 16 must be increased by a factor of 2 to obtain the 
direct- current ampere turns needed for use with a chopper. 

The ineffectiveness of the direct current to change conditions 
is further demonstrated by comparison of the voltage-wave shapes 
obtained with and without the superimposed direct current . These 
curves, which are shown in figure 17, are oscillograms from the 
search coil around the center of the blade. One set of records is 
for low values of ampere turns and the other set is for high values. 

A voltage survey along the blade (fig. 10) with this search 
coil indicated voltage variations of 14 percent at the bottom to 
-25 percent at the top compared with the value at the blade center. 
The value at the blade center used for data, was, however, approxi- 
mately equal to the average value along the blade. The inequality 
of the flux along the blade is believed to be due to the fringing 
flux from the off-center exciting coil (fig. 10). 

The large magnetomotive-force drop across the blade probably 
explains the ability of the coil fringing flux to affect the volt- 
age along the blade. The magnetomotive-force data are obtained 
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from columns 7, 8, and 9 of table I, in which the air-gap maximum 
ampere turns, the blade maximum ampere turns, and the total root- 
mean- square ampere turns, respectively, are given. These data show 
that as the blade- surf ace field intensities increase the air-gap 
ampere turns N a I max a decrease from approximately 30 percent to 

about 10 percent of the blade ampere turns Below 

6 oersteds, of course, the percentage of air-gap ampere turns is 
still larger. The large drop across the blade comes from operation 
well past the knee of the magnetization curve at the blade surface. 
Inspection of figures 4, 13, and 14 shows the knee of the magnetiza- 
tion curves to be in the neighborhood of 6 oersteds or lower, which 
is equivalent to the very low value of approximately 50 ampere turns 
on the curves of figure 16. It is recognized that this large 
magnetomotive-force drop along the blade in comparison with air 
gaps of practical sizes may make a difficult chopper- design problem) 
that is, it may be difficult to obtain sufficient chopping action at 
relatively high frequencies. 

Two curves for Armco Iron are shown in figure 16, the one pre- 
viously discussed for a total air gap of 0.04 inch and the other 
for an air gap of 0.06 inch. These two curves show that the fring- 
ing flux around the air gap is very important and, because of the 
fringing, the air-gap spacing between a tooth and a blade is not 
a critical factor in determining the maximum flux or the total 
heat generated in the blade. For tnis reason, relatively large 
clearances between blade and tooth can be accommodated if dictated 
by compressor-design criterions; however, fringing flux may also con- 
tribute to chopper- design difficulties by limiting reluctance 
change in the air gap. 

Calculated data also indicate that the flux that fringes 
around the air gap is the major portion of the flux. The calcu- 
lated permeance (table I, column 5) including the fringing space 
is approximately 9.00 centimeters) whereas the approximate per- 
meance of the air gap alone, depending upon the depth of penetra- 
tion, is only 0.5 to 1.0 centimeter. The depth of penetration is 
given in column 6 of table I and varies from about 0.005 centi- 
meter at a blade- surface field intensity E TT1nv 2 of 6 oersteds to 

about 0.015 centimeter at 50 oersteds for Armco Iron and SAE 1020 
steel. A direct result of the lower conductivity of Hipemik is 
that its depth of penetration is equal to twice that of the other 
materials. 

Equation (20a) will now be examined for the purpose of dis- 
cussing the various parameters involved. At the field intensities 


24 


NACA EM E9E06 


of interest, those that will produce 10 or 11 watts per square inch, 
the first quantity under the radical of the equation is less than 
5 percent of the radical value and can therefore be neglected. The 
equation reduces to 


P. = 


1.25/5 


sin(45°+P 2 -6 2 )10- 4 V2 cos^-k^) p 2 f 2 (H^ 2 ) 2 (B^ x ) 


lii 

)] P 2 f 2 


3 

i2 


1 

,2 


(20c) 


Also, at the field intensities of interest the values of Og and 
|3 2 -5 2 change very little, which can be seen from columns 3 and 4 
of table I for values of 2 of 30 and 50. A noteworthy fact 

is that the angle (45°+3 2 -6 2 ), the sine of which is the power 

factor, changes only from 45° to about 55° over the entire range 
of field intensity H^y g of 0 to 50. The sine (45°+0 2 -8 2 ) and 

cos (45°+a^ ) are therefore grouped in the brackets of equa- 
tion (20c) with the other constants. Equation (20c) therefore gives 
very simply the factors that determine the ability of a material to 
produce heat from eddy currents. For a given value of field inten- 
sity g or magnetizing ampere turns, the eddy current gener- 

ated per unit surface is proportional to the square roots of the 
resistivity of the material, the saturation flux density, and the 
frequency. For a given material and frequency, the heat generated 
is proportional to the 3/2 power of the magnetizing ampere turns. 
Provided the blade is sufficiently thick to make the flux at the 
center negligible, equation (20c) shows that the heating is inde- 
pendent of blade thickness and is a function of surface area only. 
Because cooling is also a function of surface area, however, the 
important criterion is the heat generated per unit surface rather 
than the total heat, and in this report equations and results are 
therefore expressed in values per unit surface. 

Because for a given installation the highest usable frequency 
is fixed by the geometry of the compressor, for a required amount 
of heat the variables will be the saturation flux density B rrin y -j_, 

the ampere turns NI, and the resistivity of the material p, or 

constant 
/^/ p ®max 1 


"W 2 = fc* 1 ) = 
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This equation shows that a blade material having high resistivity 
and high saturation flux density is desirable in the interest of 
keeping the impressed ampere turns at a minimum, because the number - 
of ampere turns required varies inversely as the cube roots of the 
resistivity and the saturation flux density. 

The advantage of a high-resistivity material is experimentally 
shown by the curve for Hipernik in figure 16. The resistivity of 
Hipernik, which is about four times the resistivity of Armco Iron, 
is the determining factor in reducing the required impressed ampere 
turns for Hipernik relative to Armco Iron for a given heat density. 
As previously indicated, Hipernik has higher flux penetration 
because of its greater resistivity relative to the other two mate- 
rials; the greater depth of eddy currents is the mechanism explain- 
ing the larger heat density. The use of a higher-resistivity 
material should help the solution of the chopper problem previously 
mentioned. 

From the data of figure 16 and for other assumed conditions, 
the direct- current power required for producing the magnetizing 
current can be roughly estimated. For one such theoretical instal- 
lation, the magnetizing power required per square inch of blade 
surface to produce a heat value of 10 watts per square inch was 
0.1 watt, or approximately 1 percent. Calculations were then made 
to determine whether the power required for ice protection would 
exceed that available from such a system of eddy-current heating. 


Calculation of Heat Eequired 

Inasmuch as the proposed system of eddy-current heating of the 
axial-flow compressor-inlet guide vanes is dependent upon a revolv- 
ing rotor, rotor speeds had to be specified for the computation of 
the heat required. No icing tolerances were allowed. Because the 
severity of icing may prove to be a function of altitude and air 
consumption, the heat requirements for several engine conditions 
and altitudes were also studied. 

The rate of heat dissipation was computed for the conditions 
subsequently enumerated. The analysis employed was taken from 
references 9 to 11 and is presented in appendix B for convenience. 

The blade-surface temperature was assumed to be kept at 35° F and 
the temperature of the compressor-inlet air was assumed to be 0° F. 
The aircraft was assumed to be flying at a true airspeed of 275 miles 
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per hour through air having a liquid-water content of 1 gram per 
cubic meter. Heat requirements were then calculated for an axial- 
flow compressor of the 4000 -wound- thrust class for the following 
conditions : 


Altitude 

(ft) 

Engine speed 
(rpm) 

Air flow 
(lb/sec) 

0 

7600 

75 


7000 

65 


6000 

50 

5,000 

7600 

65 


7000 

50 


6000 

45 

15,000 

7600 

50 


7000 

45 


6000 

35 


The rate of heat dissipation was the greatest for the highest 
speed condition at sea level (fig. 18(a)). The average heat dissipa- 
tion for this condition was approximately 5400 Btu per hour per 
square foot, which is equivalent to a power of approximately 
11 watts per square inch of vane surface. 


Shaft Power Eequired 

The power required for protecting the inlet guide vanes under 
the design conditions is approximately 11.6 horsepower or about 
8.7 kilowatts. At sea level and at rated engine speed of 7600 rpm, 
the compressor absorbs about 5680 horsepower. Thus the power 
required for protection against ice accretion would be only about 
0.20 percent of the turbine load. In addition to the power fur- 
nished by the turbine, about 1 percent of the 11.6 horsepower 
required for protection will necessarily be supplied for excita- 
tion of the coils. The electric power for excitation can be 
supplied by batteries and generator. 


Apparent Advantages and Disadvantages 

The apparent advantages of an ice-protection system utilizing 
eddy-current heating lie in the fact that heat is produced where 
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needed with little time lag. Most of the heat is used to warm the 
surface subject to icing and not the induction air as in the 
exhaust bleedback systems . The loss in net thrust would therefore 
be small because the increase in the temperature of the incoming 
mass of induction air due to surface heating would be slight. 

If ice accumulations occur before heat is applied, de-icing 
can be accomplished because eddy-current heating will melt the 
under surface of the ice coating and permit washing-back of the 
accumulated ice by the induction air. This washing -back may prove 
to be a disadvantage, however, if the pieces broken off are allowed 
to become large enough to damage the compressor blading. The 
application of some heat before the ice is allowed to accumulate 
may therefore be necessary. 

Because no special ducting is required for eddy-current heat- 
ing, the design may be kept within the limits of the maximum 
diameter of the compressor. Eddy-current heating may have to be 
supplemented with some other system of heating, however, to pro- 
tect the lips and the surface of the inlet cowling. Corrosive 
exhaust gases need not be fed back into the compressor, however, 
because they can be contained within the walls of the inlet. 

The modifications with eddy-current heating appear feasible. 
The system of protection by means of eddy-current heating can be 
applied back through the compressor for as many stages as is deemed 
necessary. The heat requirement for successive stages would dimin- 
ish with stage number. Application of this principle, however, 
becomes more difficult with successive stages because of limited 
clearances and complexity of electric circuits. 

With eddy-current heating, the use of intensity- controlled 
heating would be possible because the controls would be simple. 

The disadvantages of the application of eddy-current heating 
are the necessity for including magnetic material in the flux cir- 
cuit and the corresponding weight increase of these materials and 
of the coil and the chopper. 


SUMMARY OF RESULTS 

From an investigation of eddy-current heating for icing pro- 
tection of axial-flow-compressor blades, the following results 
were obtained ; 
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1. Heating of axial-flow-compressor blades by eddy currents 
appeared to be applicable for protection against icing. 

2. The surface average heat requirement for the inlet guide 
vanes in conditions of icing of a jet engine of the 4000-pound- 
thrust class was calculated to be maximum at rated-speed conditions 
and equal to 5400 Btu per hour per square foot or approximately 

11 watts per square inch. Creation of heat densities of this 
magnitude for ice protection of the inlet guide vanes of a com- 
pressor appeared possible. 

3. The shaft-power requirement for the generation of eddy 
currents to produce the heat needed for protection of the inlet 
guide vanes is small, and for the proposed design is only 0.20 
percent of the turbine load. Compared to this power, the elec- 
trical power required for the magnetizing field coils, however, 
is only of the order of 1 percent of the shaft power required for 
generation of the eddy currents or 0.002 percent of the turbine 
load. 


4. In an alternating- flux eddy-current system at heat densi- 
ties necessary for ice protection, the eddy-current power generated 
is proportional to the square roots of frequency, material resis- 
tivity, and saturation flux density, and to the 3/2 power of the 
maximum field intensity at the vane surface. The blade material 
chosen should therefore have the highest resistivity and the 
greatest saturation flux density. 

5. In a pulsating-unidirectional flux eddy-current system in 
which the flux varies from zero to a maximum value and which can 
be simulated by superimposing a direct-current flux on a sine-wave 
alternating- current flux, the heat generated does not change appre- 
ciably with application of the direct- current flux. Heat calcula- 
tions can therefore be made on the basis of the alternating- current 
component. In the production of the pulsating flux by means of a 
chopper and the application of a direct- current magnetomotive force, 
however, the required direct- current ampere turns will be greater 
than the root-mean- square ampere turns of the alternating- current 
component by a factor of 2 /\J2. 

6. The vane surfaces will probably operate well past the knee 
of the magnetization curve of the blade material, and the 
magnetomotive-force drop along the blade will be great in com- 
parison with air gaps of practical sizes, which may lead to diffi- 
culty in obtaining sufficient chopping action at high frequencies. 
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7. Provided that the blade is sufficiently thick to insure a 
negligible flux at the center, the heat generated per unit surface 
is independent of vane thickness or the total heat generated is a 
function of only the vane surface. 

8. The disadvantages anticipated in the application of eddy- 
current heating would be the necessity for including magnetic 
material in the flux circuit and the corresponding weight increase 
of these materials and of the coil and the chopper. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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APPENDIX A 


DERIVATION OF ANGELES & 2 AND 3 2 
From figure 6, 

d6 = B max 1 62 sin ( 45 °+3-S) 


$> 


max 


(Al) 


or, using equation (10), 


d6 = 


Brnny i dZ sin (45°+3-6) 


-^max 1 
a/2 c 


Bmax 1 d * cos (45+3-6) 


sin(45°+3-6) ^ c dl cos (45°+3-5) 

cos (45°+3-6) 


(A2) 


1 + / a/2 c dZ cos (45°+3-6) 

Integrating in the limits 0 to 5 2 and 0 to Z 2 , 

6 2 = tan (45°+02)log e jl + a/2 c cos (45°+02) Zgj (A3) 

where Og is the mean value of Pg -5 2 for the S iven value of *2> 
which gives from equation (15a) 


& 2 " 


tan (45°+02) 


log c 


1+2 cos (45°+o 2 )( 


,x 2 


v^raax 1 


- 1 


(A4) 


The derivation of 3 2 follows a similar procedure. In fig- 
ure 6, the increment of current in each new sheet of material of 
thickness dZ is equal to the product of induced voltage and the 
conductance of the sheet, or y E ^y dZ ; therefore. 


I = I 
max max 


i + y E dZ sin (45°+3-&) 

I J max 

= “TOT + f 7 E max dl sin (4s°+e-&) 


(A5) 
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From the vector diagram, the following relation is true: 
I max d P = 7 Emax dl cos (45°+p-8) 
or, when equations (A5) and (A6) are combined. 


(A6) 


0.4 icy E 


max 


dp = 


cos (45°+3-6) Emajj q 


dZ sin (45°+p-6) 


sin (45°+P-8) po.4 ny E 


(A7) 


1 + 


max 


^max 1 


dZ sin (45°+p-6) 


Integrating from 0 to P 2 and from 0 to 1*2’ w ^ ere is 

the average of 3 2 -6 2 as before, yields 


P 2 = cot (45 0 +02)log e 


1 + 


0.4 *7 Emax , o 

dZ sin(45 +p-8) 


10 


^max 1 


— aA8) 


The first step in the evaluation of the integral is to replace 
E mqy by using equations (10) and (12). 


0.4 «7 
Emax 1 


dZ sin (45°+3-6) = 


0.4 try 2 rtf 10~ 8 


10 

or 


®max 1 


J2 


+ B max 1 1 cos (45°^) 


dZ sin (45 +3-5) 


Jo 


tj2 c + 2c 2 Z cos (45°+a) 


dZ sin (45°+3-8) 


or where a is again the average value of (3-6), 


0.4 icy EL-^ r _ o + 

— H mflV 7^ dI sin (45°+3-5) = \jz cl 2 + c l 2 cos(45 0 +a2) 


(A9) 


sin (45°+a2) 
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With substitution of equation (A9), equation (A8) becomes 


3 2 = cot (45° 


-KXg ) log e il + 1^2 c + c 2 2 2 cos( 45 0 +02)J sin(45°+a2)) 

^ (A10) J 


Substitution of equation (15b) into equation (A10) gives the 
solution for 


P 2 = cot (45°+a 2 ) log e 


+ S in(45°^)^-l) 


(All) 
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APPENDIX B 


CALCULATION OF RATE OF HEAT REQUIRED TO MAINTAIN SPECIFIED 
BLADE-SURFACE TEMPERATURE IN CONDITIONS OF ICING 

An analysis of the heating requirements for thermal-ice pre- 
vention of the inlet guide vane of an axial-flow compressor is 
presented for specified conditions of icing. The feasibility of a 
design of an ice-prevention system utilizing eddy-current heating 
under the assumed icing conditions has been presented in the text 
of this report. 


Symbols 

The following symbols are used in this appendix: 

A cross-sectional area of flow passage, sq ft 
C blade chord, ft 

Cp specific heat of air at constant pressure, Btu/(lb)(°F) 

Cp w specific heat of water, Btu/(lb)(°F) 

D diameter of equivalent leading-edge cylinder, ft 
d average distance between adjacent blades, ft 
e vapor pressure, lb/sq in. 

e g vapor pressure at saturation at temperature of surface, 

lb/sq in. 

e x local vapor pressure, determined by equation (B9), lb/sq in. 

g acceleration due to gravity, 32.2 ft/sec^ 

H total heat dissipation, Btu/(hr)(sq ft) 

H-l heat dissipation due to convection, Btu/(hr)(sq ft) 

H2 heat dissipation due to evaporation, Btu/(hr)(sq ft) 

H 3 heat dissipation due to heating of impinging water on surface 

of blade, Btu/(hr)(sq ft) 
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h unit thermal convective conductance between vane surface and 
air, Btu/(hr)(sq ft)(°F) 

J mechanical equivalent of heat, ft-lb/Btu 

L latent heat of vaporization of water, Btu/lb 

Z length of blade, ft 

M water impingement per unit area, lb/(hr)(sq ft) 

Mg water impingement per unit area leading edge, lb/(hr)(sq ft) 

M s water impingement per unit area blade surface, lb/(hr)(sq ft) 

m liquid- water content, lb/cu ft 

N number of blades 

Pr Prandtl number 

p static pressure, lb/sq in. absolute 

R gas constant, for air 53.3 ft-lb/(lb) (°R) 

s distance along chord from leading edge of blade, ft 

T av average of inlet-air temperature and blade- surface tempera- 
ture, °R 

t air temperature, °F 

t g surface temperature with no external heat conditions, °F 

t s ^ surface design temperature, °F 
V velocity, ft/sec 

W air flow, lb/sec 

x average distance between adjacent blades, ft 

a angle of attack of blade section, degrees 

7 ratio of specific heats, for air = 1.4 

p density, slugs/cu ft 
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Subscripts: 

c compressor inlet 

s surface condition 

x local condition at any chordwise position on blade surface 


Liquid Water in Atmosphere 

The basis for the design of an ice-prevention system was 
established from recent and extensive information on the severity 
of icing conditions likely to be experienced by an airplane under 
normal engine operation (reference 9). The icing condition selected 
may occur in cumulus clouds having a liquid-water content of 
1.0 gram per cubic meter and an average drop size of 20 microns in 
diameter. The inlet temperature was assumed to be 0° F. 


Impingement of Water on an Airfoil 

Little information is available regarding the impingement of 
water on airfoils as small as those found at the inlet of an axial- 
flow compressor of the type studied herein. Limited use was there- 
fore made of the existing theory involving much larger airfoils, 
such as propeller blades and wings. Figure 19 indicates the shape 
and dimensions of the inlet-guide vanes considered in this report. 
Eighty-eight such blades are at the compressor inlet. The water 
impingement was assumed to be confined to the leading edge and the 
concave surface of the vane and, because the vanes are small, their 
collection efficiency was further assumed to be 100 percent. 

The specific conditions for which the rate of heat required is 
calculated are for a compressor- inlet air temperature of 0° F and a 
vane-design surface temperature of 35° F. The temperature and the 
pressure of the air passing over the first row of vanes differ from 
the inlet-air temperature and pressure because of kinetic effects 
and are functions of air velocity and density. 

Velocity gradient through blading. - Any two adjacent blades 
of the same compressor stage, because of their shape and spacing, 
can be assumed to constitute a nozzle. The cross-sectional areas 
used for computation of the average air velocity over each section 
of blade can be determined using average values of the distance 
between adjacent blade surfaces. 
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The section flow area is 


A x - *^x 


(Bl) 


The spanwise average air velocity through a passage at any 
blade section is then 


V x p xgAx N 


(B2) 


Pressure gradient through stage. - The air pressure in the 
passage between the blades is determined assuming a reversible 
adiabatic-state change expressed in terms of the initial state and 
velocity at the blade section. 


Px = 




2gRT c 


7 

7-1 


(B3) 


The velocity V x as determined from equation (B2) 'is a func- 
tion of density p x and area A x . Inasmuch as p x is a function 
of p^, solution of equations (B2) and (B3) would involve exten- 
sive computation. In order to simplify the work, p^ was there- 
fore assumed equal to p . 

c 

Surface temperature in dry air. - Hardy (reference 10) has 
indicated that in clear air flow the temperature of the unheated 
surface for laminar flow will equal 



and for turbulent flow will equal 
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The index of Pr is dependent upon the type of flow, laminar 
or turbulent; for the conditions of icing used in this analysis a 
Prandtl number of 0.72 was assumed. 


Unit thermal conductance. - Three values of unit thermal con- 
ductance taken from reference 11 were considered in this analysis - 
those for stagnation, laminar, and turbulent flow. 


Stagnation: 


Laminar flow: 


h 

x 


= 0.194 T 


av 


0.49 





= 0.0562 T av 


0.5 



(B6) 


(B7) 


Turbulent flow: 


0.524 T, 


0.296 


av 


' V x«>c 

„0.25> 


0.8 


(B8) 


The temperature T QY is the average of inlet-air temperature and 
blade-surface temperature. 


The location of the point of transition from laminar to tur- 
bulent flow is important because of its effect on the heat- transfer 
coefficient. Because there is no reliable method of determining 
this point with the presence of water on the surface, however, a 
location of 10-percent chord was selected. 


Local vapor pressures. - The assumption is made that the air 
is initially of high humidity and, because of a rapid change in 
pressure as the air passes over the vanes, there is little or no 
change in state except that the air becomes supersaturated. The 
local vapor pressure will then be 


e 


x 



(B9) 
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Heat required. - The heat dissipation from a surface in con- 
ditions of icing may he divided into its four components: 

(l) Dissipation of heat due to convection from the surface to 


(2) Dissipation of heat by evaporation of water from surface 

because of difference in vapor pressure 

(3) Dissipation of heat as a result of raising temperature 

of impinging water 

(4) Dissipation of heat due to radiation from relatively warmer 


For this analysis, the convex surface of the blade was assumed 
to dissipate heat only by convection and evaporation, because very 
little, if any, water is assumed to impinge upon the convex surface 
of the blade. Heat losses due to radiation were small and were 
neglected. 

Dissipation of heat due to convection and evaporation. - The 
heat dissipated as a result of convection and evaporation was 
computed by the following simplified equation, as derived from 
reference 10 : 


Water impingement. - The water impingement per unit area of 
leading edge and vane surface was computed assuming a collection 
efficiency of 100 percent. The areas affected were evaluated from 
the geometry of the guide vane with straight-line impingement 
assumed on all exposed surfaces. 

The average water impingement per unit area of leading edge 
was calculated as follows: 


the air 


surface 



(BIO) 


M_ = - Y r m (3600) 


(Bll) 


The water impingement per unit area of vane surface was calculated 
by the following equation, which represents the concave surface as 
a plane surface parallel to the chord line: 
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M g = 3600 m V c sin a (B12) 


Dissipation of heat as result of raising temperature of 
impinging water. - At the leading edge and for the concave "blade 
surface, the rate of heat dissipation was computed with the assumed 
collection efficiency. The heat required to raise the temperature 
of the impinging water to that of the blade was therefore equal to 

% = MCpw (t s d - t c ) (B13) 


neglecting the small kinetic temperature rise in the inpinging 
water . 


Total heat dissipation. - The total heat dissipation for the 
leading edge and the concave blade surface is the summation of the 
three types of heat loss 


H = (VV + =3 


(BU) 


whereas the total heat dissipation of the convex blade surface was 
assumed to be 


H = H-^+Hg 


(B15) 


The average heat dissipation presented per square foot of vane 
surface (fig. 18(a)) is the integrated average of the dissipation 
rates for both sides of the vane. The heat dissipation due to 
evaporation, Hg in equation (B15), is assumed to apply because some 
moisture will probably run back on the convex face. Also, the heat 
requirement for the convex face must be nearly as high as that for the 
concave face to prevent the development of regions of subfreezing tem- 
peratures on the blades. The increase in heating requirement that 
results from this conservative assumption should help to safeguard 
the leading edge by action of conduction in the blade metal. 
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TABLE I - SUMMARY OP PERTINENT VALUES 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ii 

12 

13 

14 

15 

H max 2 
(oer- 
steds) 

H raax 1 
(oer- 
steds ) 

°2 

(deg) 

P 2 " 6 2 

(deg) 

p T 

(cm) 

Zqt 

(cm) 

^a^max a 

^b^max b 

MI 

^2 

(volts/ 

cm) 

Measured 

E 2 

(volts/ 

cm) 

V 

(watts/ 
sq in. ) 

Measured 

P 1 

(watts/ 
sq In. ) 

n 

c 

(1/ cm) 






SAE 1020 steel. 

measured y 3 

0.594 xlO 5 mho/ c 

■m 




6 

4.31 

0.60 

1.30 

8.60 

0.0049 

16.8 

56.1 

51.6 

0.0140 

0.0127 

0.22 

0.25 

2665 

195.5 

30 

5.63 

4.25 

8.62 

9.01 

.0123 

37.2 

280.5 

224.7 

.0379 

.0400 

3.33 

3.50 

2665 

195.5 

50 

5.93 

5.26 

9.82 

9.24 

.0163 

46.7 

468.0 

364.0 

.0502 

.0535 

7.44 

7.60 

2665 

195.5 




Annco Magnetic Ingot Iron, measured y 

* 0.873 xlO 5 mho/ cm 




6 

2.83 

1.65 

3.90 

8.56 

0.0043 

14.7 

56.1 

50.0 

0.0128 

0.0115 

0.21 

0.21 

4590 

311 

30 

3.49 

5.30 

9.88 

8.90 

.0104 

31.1 

280.5 

220.2 

.0323 

.0350 

2.87 

2.80 

4590 

311 

50 

3.60 

6.05 

10.10 

9.10 

.0138 

40.0 

468.0 

359.0 

.0426 

.0444 

6.34 

6.70 

4590 

311 






Hipernlk, measured ^ 

• = 0.: 

248 x 10 5 mho/ cm 





6 

0.256 

6.68 

10.15 

8.66 

0.0104 

23.7 

56.1 

56.4 

0.0241 

0.0240 

0.43 

0.40 

50,000 

546 

20 

.280 

7.55 

10.25 

9.13 

.0190 

42.6 

187.0 

162.3 

.0458 

.0510 

2.74 

2.50 

50,000 

546 

40 

.300 

7.88 

10.30 

9.54 

.0264 

59.3 

374.0 

306.5 

.0667 

.0780 

7.96 

8.40 

50,000 

546 


NACA RM E9E06 




7 ^-/ 4 / 2 . . . 1125 


NACA RM E9E06 


43 


r 1 



U A 



\\\\\ Vv 





r 

l 

I 



Figure 2. - Schematic diagram illustrating principle of eddy-current heating 
applied to inlet guide vanes of axial-flow compressor. 
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Figure 3. - Section of plate. 
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Figure 4, - Normal magnetization curve for SAE 1020 steel annealed at 1500° F and furnace -cooled* 

(Data from reference 5 # ) 
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-E mflT 3 , voltage consumed 
by flux s 


Exoax s t voltage 
Induced by 
flux a 



Figure 5. - Vector diagram of phase relation of current, flux, 
and voltage at constant permeability. 
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Figure 6. - Vector diagram of phase relation of current, flux, and 
voltage under conditions of saturation. 
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Figure 8. - Variation of a 2 with 2 /Hma x 1 
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Figure 9. - Variation of P2~®2 with H m9 Y 2/^max 1 
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Figure 10. - Equipment utilized for determining eddy-current heat generated 

in blade. 
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Figure 11. - Circuit schematic diagram for alternating-current 

plus direct-current experiment. 
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Figure 13. - Normal magnetization curve for Armco Magnetic Ingot Iron. 

(Data from reference 7. ) 
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Magnetic-flux density, gausses 
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Figure 14. - Magnetization curve for Hipernik, 
(Data from reference 8. ) 
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Total 

Calculated 

Measured 


air gap 

cl • C • 

a.e# 


(in. ) 


o 

SAE 1020 steel 

0.04 

X 

V 

Armco Magnetic 

.04 


Ingot Iron 


A, 

A 

Hipernik 

.04 



0 100 200 300 400 


Effective ampere turns 

Figure 15. - Variation of voltage of self-induction per centimeter of blade width with 
ampere turns. Vane, 4>| by 1 by g Inch (Hipernik > 4g by 1 by 0.094 in.); frequency, 
6100 cycles per second. 
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Figure 16 


# - Ampere turns required for eddy-current heating# Vane, 4| by 1 by 
(Hipernik by 1 by 0.094 in.); frequency, 6100 cycles per second. 
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(a) 99 root-mean- square ampere turns, (b) 89 root-mean- square ampere turns; 139 direct- 

current ampere turns. 
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Figure 17. - Search-coil voltage oscillograms for comparing wave shapes with and without superimposed direct- current ampere 

turns. Armco Magnetic Ingot Iron; air gap, 0.04 inch. 
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(c) 383 root-mean- square ampere turns. 


(d) 281 root -mean- square ampere turns; 414 direct- 
current ampere turns. 



C- 23406 
5-6- 49 


Figure 17. - Concluded. Search-coil voltage oscillograms for comparing wave shapes with and without superimposed direct- 

current ampere turns. Armco Magnetic Ingot Iron; air gap, 0.04 inch. 
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Heat dissipation, H, Btu/(hr)(sq ft) 
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(c) Engine speed, 6000 rpm. 

Figure 18. - Calculated unit heat dissipation along inlet-guide-vane chord in icing 
condition. Blade-surface temperature, 35° F; inlet-air temperature, 0° F; liquid- 
water content, 1 gram per cubic meter. 
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Figure 19. - Physical dimensions of inlet guide vanes of typical axial-flow 

compressor. 
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